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ABSTRACT 

Relations for the asperity-contact time fraction during elastohydro- 
dynamic lubrication of a ball bearing are presented. The analysis is 
based on a two-dimensional random surface model, and actual profile 
traces of the bearing surfaces are used as statistical sample records. 
The results of the analysis show that transition from 90 percent contact 
to 1 percent contact occurs within a dimensionless film thickness range 
of approximately four to five. This thickness ratio is several times lar- 
ger than reported in the literature where one-dimensional random sur- 
face models were used. It is shown that low pass filtering of the statis- 
tical records will bring agreement between the present results and those 
in the literature. 

INTRODUCTION 

Elastohydrodynamic lubrication (EHD) is the term used to describe 
part of the technology concerning lubrication of concentrated mechanical 
contacts. In essence, EHD lubricant film formation depends on the 
coupled effects of physical changes in the lubricant, which are caused by 
high pressures in the Hertzian contact area, and elastic changes in the 
shape of the Hertzian contact area, which affect the pressure distribu- 
tion . The high pressures in the EHD contact area act to squeeze out the 
lubricant. However, the lubricant becomes thicker (more viscous) with 
increasing pressure and resists being squeezed out. The net result is 


the formation ol a thin lubricant film that is beneficial in preventing 
seizure and rapid wear of the contacting parts (refs. 1 to 4) . For many 
applications the EHD film thickness is the same order of magnitude as 
the surface rms roughness. Experimental measurement of the film 
thickness is very difficult because films are so thin. Various methods 
that have been used are optical (interferometry), X-ray, and electrical 
capacitance and conductance techniques (ref. 5). Of the aforementioned 
measurement methods, the capacitance and conductance methods are 
most suited to measurement of film thickness in full-scale bearings 
(refs. 6 and 7). The conductance method of measurement depends on 
having a known relation between film thickness and contact time frac- 
tion, The contact time fraction is directly related to the normalized av- 
erage voltage observed when a low voltage is applied across the lubri- 
cant film. 

In 1964, Tallian and his coworkers (ref. 7) formulated a statistical 
model of bearing surface roughness and used the model to infer EHD 
film thicknesses, based on electrical conductance measurements. Their 
results were applicable to the regime of "partial EHD contact," where 
the load is shared by the EHD film and the high points or asperities of 
the metal surfaces that momentarily interrupt the lubricant film (ref. 8), 

By the early 1970’ s it was generally accepted that asperity contact 
must be viewed as a random process (refs, 8 to 14). Most researchers 
used stylus traces of the surface to obtain profile statistics for the ran- 
dom process models. In 1971, Nayak (ref. 14) explained how Longuet- 
Higgens' theoi’y of ocean surfaces (refs. 15 to 17) could be used to model 
rough surfaces as two-dimensional, isotropic, Gaussian random pro- 
cesses. He showed that significant differences exist between surface 
statistics and profile statistics and that a naive analysis assuming that 
profile statistics may be directly used is erroneous (refs. 14 and 18), 
Sidik has extended the theory of Nayak to obtain a model for asperity- 
contact time fraction as a function of film thickness in partial EHD con- 
tact lubrication (ref, 19). In reference 20, the theory Is generalized to 
nonisotropic Gaussian surfaces. In reference 21 the computational 
methods are developed. 

The objective of this paper is to apply the relevant results of two- 
dimensional random surface analysis to obtain a relation between 
asperity-contact time fraction and average EHD film thickness for a 
typical ball bearing and to compare the results with those expressed in 
references f> and 11. Also, recent work by Sayles and Thomas (ref. 22) 
is used to correlate the results. 
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Asperity- Contact Model 


The bail bearing for which this analysis was performed has a 20- 
mllllmeter bore and three 7. 15-mlllimeter- (9/32-ln,-) diameter balls. 
Three different thrust loads were considered in the analysis of contact 
time fraction. Table I gives the calculated Hertzian stresses and con- 
tact ellipse dimensions corresponding to the different loads. 

Under loaded conditions, assume that the ball and race surfaces are 
two-dimensional ergodic Gaussian processes and that within the Hertzian 
contact zone the mean planes are parallel and separated by a lubricant 
film of thickness h. A cross section of a single ball- race contact is 
presented in figure 1. Coordinate x is in the direction of rolling. The 
ball surface is denoted by z^^^.y) and the race surface by Zj.(x,y). The 
two processes and z^ are independent, with mean levels = 0 

and Pj. = 0, correlation functions and R^, and variances and 

(7^, respectively. 

The composite process z = Zl + z_ is also an ergodic Gaussian 

or o 

process with mean zero and correlation function R and variance , 
where 


V = ogR^(r^, Tj,) H- Ty) (1) 

With this notation, then (as shown in fig. 2) any metallic contact occur- 
rence is represented by the composite surface rising above the level h. 

An approximation to the time fraction during which there is metallic 
contact anj^here within the Hertzian zone will now be derived. Consider 
the process z above the x-y plane. At a level h above the reference 
plane, pass a cutting plane that will occasionally intersect z. The sets 
of points in the reference plane where z(x, y) 2 : h are called excursion 
sets. Such excursion sets are represented as the crossed areas in fig- 
ure 3, Superimposed upon this plane is an elliptical region that repre- 
sents the Hertzian contact area. For constant rolling velocity, this el- 
liptical region moves to the right at a constant velocity v through the 
region bounded by the parallel dashed lines y^ and y^. At the termi- 
nation of a test period of time T the Hertzian area is at the elliptical 
region at the right in figure 3. 

If it is assumed that A = h/a is large there will be few excursions 
of z above h and, hence, few metallic contacts. The contacts will be 
small in area, and the probability of two or more contacts in a small 
area is negligible. The contact occurrence is as follows: The dashed 
ellipse on the left in figure 3 represents the location of the Hertzian area 
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when the contact Is first made. The dashed ellipse on the right repre- 
sents the location of the Hertzian area when the contact is broken. The 
two points and P(Xg,yg) denote coordinates of the make- 

contact and break-contact occurrences. The distance of P(Xg»yg) to the 
centerline of the ellipse on the right is termed Lg. Thus, contact exists 
for a total distance L that is composed of three parts . Two of these 
parts are Lj^ and LgS the third is termed X and is the distance 
Xj^ - Xg. It is assumed that X is negligible with respect to and 
Lg and the excursions are uniformly distributed with respect to the 
y-axis. Also Lj^ and Lg are approximately equal. As a result, the 
statistical expectation of the contact distance L is 

E{L} ^2E{Lj^^} (2) 

where Trf/2 is the average length of the Hertzian ellipse. The expected 
total contact time EfT*} can be approximated by the product of the av- 
erage number of excursions and the average time of contact for each. 
Equivalently, E{T*} can be approximated by the expected number of ex- 
cursions per unit area E{\} times the area rolled over, times the av- 
erage time of contact for each. If T is large, the rolled-over area is 
approximately vTw, so that 

E{T+} = E{x}(yTw)^^ (3) 

The expected contact time fraction is obtained as 

E{T^} - iilEi = iiiil E{x} (4) 

T 2 

Thus, E{T^,} as a function of A can be calculated directly from a com- 
putation u*’ E{x} as a function of A. The determination of E{\'} is 
presented next. 

Following the method of Nayak (ref. 14), let the following variables 
^|(i - 1,6) be defined; 
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4l(x,y) - -^z(x,y) 
8x^ 


U - 4o(x.y) - -^z(x,y) 
^ Ox 8y 


= ^3(x.y) = -^z(x.y) 

9y^ 


^4 = 44(x,y) - z(x,y) 


4r, = ^(x,y) = — z(x,y) 
^ Ox 


= — z(^»y) 

^ Oy 




( 5 ) 


The expected number of summits of height within a unit area Is 
given by the triple integral 


f«4) 


|y|l/2|sil/2p^j3 





( 6 ) 


where the region of integration V is defined by 


^ 3^0 






( 7 ) 


and and S by the moments of the power spectral density (see 
appendix) 
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Equation (6) is transformed to cylindrical coordinate:' followed by a ro 
tation. The transformation equations are 


t - r - p cos ffi 

t g-gi!Lg 

t = r P cos 0 

^3 r 

\J2 


t - A ml/2 

^4 ~ ^“OO 


ai) 


It now becomes evident that the region V describes a semi-infinite cone* 
as shown by the limits of integration on the transformed equation . 
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Equation (12) was evaluated numerically on the digital computer. The 
expected number of summits per unit area is given by the in- 

tegral 


sum 


= / f(A)d(A) 


(14) 


The probability density for summit heights is given by the ratio 


p*(A)=iI/iL (15) 

^sum 

The expected number of excursions above level A per unit area is ap- 
proximated by the product of peaks per unit area and the proportion of 
such peaks that exceed the level A . 


E{X(A)} 



P*(OcU 



( 16 ) 


This approximation is valid only in the limiting sense as A 

Adler and Hasofer (ref. 23) also provide an approximation for the 
uperossings of a process z(x,y) over the level A. In terms of the pres- 
ent notation, their results are given by the relation 
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E{x(A)} « 


AlSl^/2 (: 


(2„)3/V 


(17) 


It is interesting to note that this expression does not involve any fourth- 
order moments. 

The derivation to this point provides contact time fraction as a func- 
tion of h for a single Hertzian contact. Next, the results are applied to 
a ball bearing with three balls for three different loads. It is assumed 
that at each of the six ball-race contacts the mean film thickness is the 
same and that each contact is statistically independent of the others. 
Because of the geometry of the balls and the race, however, the nomi- 
nal Hertzian areas at tlie inner and outer race contacts are different. 
Table I presents the calculated conditions at these contacts for three dif- 
ferent loads. From equation (4) it is evident that l.{T } is simply the 
area of the Hertzian contact times E{x} . 

Let T^ and T^ denote the expected contact time fractions 
at the inner and outer races for a single ball. Thus the probability of 
no contact on a single ball is 1 - .InTc.ouf From the independence 
assumption, the probability of no contact on any of the three balls is the 
quantity 


T, 


1 - T 


c,ov 


= (1 - T 


c, in ^c, out 


(18) 


Surface Measurements 

The first step in determining the expected number of excursions per 
unit area is to obtain and analyze surface profile traces from the bear- 
ing surfaces. By using these profile traces, the important surface pa- 
rameters are computed as outlined in the previous section. As ex- 
plained in the appendix, in order to characterize the surface statistics, 
one must first obtain profile traces in at least three different direc- 
tions. The surfaces of the ball and race specimens were traced to obtain 
records of the surface profile. Several traces in different directions on 
the ball surface showed that it was an isotropic surface; the race surface 
w'as nonisotropic. 

It was not possible to obtain a sample record on the race surface in 
any direction other than the rolling direction and the cross groove di- 
rection. In order to obtain the necessary additional traces, a flat speci- 
men was prepared by material and finishing methods identical to those 
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used in making the bearing race. The flat specimen was approximately 
2.5 by 5,0 centimeters (1 by 2 in.) with a 0.13- to 0 . 25-micrometer 
(5- to 10-/iin.) CI^ surface finish. Similarly, a surrogate ball, 14.3 
millimeters (9/16 in.) in diameter with a 0.03- to 0,05-micrometer 
(1- to 2-nln.) CLA surface finish, was used to obtain sample ball rec- 
ords of sufficient length. 

Figures 4 and 5 show the microtopography of the ball and the flat 
specimen. Traces to be analyzed were recorded for three different di- 
rections on the ball and six different directions on the flat specimen. 
Nominally, the traces were taken 45° apart for the ball and 18° apart 
for the flat specimen. The first trace on the flat specimen was taken 
in the direction of the lay of the surface finish. This Is the rolling di- 
rection for the ball in the race. Five more traces were taken, with the 
last trace being at 90° to the lay of the surface finish. The traces were 
recorded in analog form as an FM signal on magnetic tape. This tape 
was sampled at equal intervals and written in digital format on another 
tape. The sampling intervals and the number of sampled points, along 
with the total sampled length, are presented in table n. 

Each of the digitized profile traces was processed by a computer 
program that performed the following steps; 

Plotted input . - All the plots were examined for faithful reproduc- 
tion of the original record and were visually found to be Identical. 

Moving-average trend removal . - The raw data required detrending 
for two reasons. First, the stylus head does not follow a path parallel 
to the mean line ofth j profile when tracing. This causes a linear trend. 
The other reason foi' detrending was to anticipate the detrending that 
occurs naturally in the lubrication process. Therefore, it was decided 
to remove trends with wavelengths longer than the Hertzian contact. 

The moving-average trend remover, which is essentially a high-pass 
filter, was effective in removing these trends. 

The number of points ir. the moving average is a function of the sam- 
ple interval and the load since the contact ellipse dimensions change with 
load. The number of points for each load and each profile are given in 
table HI . 

Normal probability plot . - A simple graphical test for normality is 
to plot the sample cumulative distribution function on Gaussian probabil- 
ity paper. The plots indicate general agreement with Gaussian distri- 
bution (fig. 6). 

Estimated spectral moments . - Spectral moments are estimated by 
the variance of the derivatives where the derivatives are approximated 
by the finite difference method. Table IV presents the estimated mom- 
ents for the three load conditions. 
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Adding ball nwments to flat profile moments . - According to the de- 
finitional equatloi. for spectral moments (eq. (A8)), the spectral mom- 
ents of a prof Hi; of the composite surface are simply the sums of the 
appropriate individual surface moments. 

Estimating spectral moments by least squares . - The following re- 
lations between the two-dimensional moments m^j and the profile mom- 
ents ^ are obUiined from equation (A9), 

2 2 
"^2 0" *”20 ® 0 sin 0 + mp 2 sin 6 (20) 

4 3 2 2 

*”4 0 ~ *”40 ® ^*”31 0 sin 0 + 6 m 22 cos 0 sin 0 

3 4 

+ 4mj^g cos 0 sin 0 + m^^ sin 0 (21) 


Equation (19) implies that the best estimator for m^^ is simply the av' 

-2 -2 

erage of the m^ The results for m^Q are 5.97x10 , 6.41x10 , 

.9 

and 5.90x10 square micrometer, respectively, for loads of 90 , 445, 

and 3100 newtons (20, 100, and 700 Ibf). 

Equation (20) provides one equation for each 0, or a total of seven 

equations in three unknowns. For each load the estimates for m 2 Qi 

mj^|, and m^g sre given here in a matrix format, emphasizing that 

this represents the covariance matrix of <p and (p in the spectral 

X y 

density function. The matrices are nearly diagonal, which means that 
the distributions of and <Py are effectively uncorrelated. For the 
90 -newton (20 -Ibf) load, 


*"20 

mil 


2.39x10'^ 
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mil 

*"02 


3.4x10"^ 

5.73x10“^ 







For the 445-newton (100-lbf) load, 


•• 





*”20 
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1.31x10"'* 
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1.31x10"'^ 

5.58X10“^ 


( 22 ) 


( 23 ) 
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For the :?iOO-nt'\vton (700-lbf) lo.id, 


m.,0 


2.17X10"'^ -4.4X10"® 

"'ll "’02 


-4.4x10"® 5.66x10“^ 


( 24 ) 


Kquation ( 21 ) provides one equation for each 0 , hence, seven equations 
in the five unknowns ni.jj, *^ 13 ’ ^ 01 ‘ f'®*' each load 

the estimates of these moments are as follows; For the 90-newton 
( 20 ”lbf) load, 
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(25) 


For the 445-newton (lOO-Ibf) load, 
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RESULTS AND DISCUSSION 

The final result of this analysis was obtained by using the moment 
estimates for the two-dimensional process to calculate the expected num- 
ber of excursions per unit area (eqs. (12) to (16)) for various film thick- 
ness ratios . These expectations were then used to obtain the no-contact 
time fraction T^ (eq. (18)). The results are presented in table V. 

Curves of no-contact time fraction against film thickness ratio are 
plotted in figure 7. The cslcu’^tions for these results required some 
30 hours of computer time. 

For comparison, values of £{x| were calculated from the spectral 
moments in table IV by using the method of Adler and Hasofer(eq. (17)). 
Tiie resulting no-contact time fraction for the 90N (20 lb) load is also 
plotted in figure 7 for comparison. The spacing between curves for the 
445N and 3100N loads was similar, but are not shown on figure 7. The 
Adler and Hasofer method gives curves that are shifted toward lower 
values of A by an amount approximately equal to 0.35. The reason for 
this difference is unclear. Both methods should be asymptotically 
equivalent, and values of A greater than 4 should be sufficiently large 
for the asymptotic results to hold. Adler and Hasofer* s "upcrossings" 
seem closer to the required "excursion sets" in the development of 
contact time fraction than the approximation based on peak height distri- 
bution. Nevertheless both arc approximations. 

Measurements of the no-contact time fraction by Talllan, et al. 

(ref. 7) on a 4-ball apparatus and by Poon and Haines (ref. 12) on a 
point contact disc machine are shown in flgure 7 for comparison. The 
measurements indicate that the no-contact time fraction begins to in- 
crease from zero at A ^ 1 and monotonically increases until there is 
100-percent film at A ^ 3.5. By comparison the results of this inves- 
tigation show this change in contact time fraction occurring for film 
jiicknesses several times larger. Also, the incremental change in 
nim thickness corresponding to the incremental change from no-contact 
to 100-perccnt contact is smaller. 

The theoretical differences in probability distributions for peak 
heights of summits on a two-dimensional surface and peak heights on 
profile traces alcmg a fixed direction on the surface may partially ac- 
count for these differences. However, the recent work of Sayles and 
Thomas (ref. 22) provides much more insight. They have shown that the 
asperity density and spectral moments are strongly dependent, on the 
surface sample interval. Their work shows that an order of magnitude 
decrease in sample interval causes approximately an order of magni- 
tude increase in the second-order moments and two orders of magnl- 
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tudf increase in the fourth-order moments. Also, as the sample Inter- 
val becomes smaller than 10 pm, the summit density seems to increase 
without approaching an asymptote. Therefore, it was decided to com- 
pare the present asperity-contact model to several models in the litera- 
ture within the framework of effects of asperity summit density and order 
of magnitude changes in the spectral moments. Basically these ar^ the 
effects of low pass filtering of the surface traces. 

The surface traces used for the work in this paper were sampled at 
0. 31 pm intervals which was sufficient to accurately reproduce the ana- 
log traces. It is not practical at this time to filter and reprocess the 
original surface traces. However, the effect of filtering can be pro- 
duced by using Adler and Hasofer’s result with the second-order mon - 
ents reduced by one or two orders of magnitude. The result of this pro- 
cedure is shown in figure 7 . The effect is to bring the T^ curve closer 
to the result generally reported in the literature. 

Curves were calculated using the method of Johnson, Greenwood, 
and Poon (ref. 11), and Kannel and Snediker (ref. 6) which is based on 

an earlier paper by Greenwood (ref. 9). Reference 11 > -umes that 

2 

there are 20,000 summits/cm and reference 6 assumci '0,000 summits/ 

2 

cm where the assumption was "used to obtain a reasonable fit with the 
experimental results" (ref, 11) which were reported by Tallian (ref, 7) 
and Poon and Haines (ref. 12). 

Results were also obtained by substituting the values of E{i^} from 

table V into Johnson, Greenwood, and Poon’s Poisson process model. In 

addition, results were obtained using the Kannel and Snediker model but 

8 2 

it is now assumed that there are ~10 summits/cm . The results were 
identical for the two methods and there was exact agreement with the 
method of Coy and Sidik near the full film conditions as shown in fig- 
ure 7. 

P’rom the comparisons presented it is clear that the main reason for 
the disparity in results is due to the effects of sample interval which is 
equivalent to low pass filtering. 

Some comments are needed regarding the usefulness of these curves 
as a means of measuring film thickness by the electrical conductance 
method. The results presented show the incremental change in film 
thickness corresponding to the incremental change from no contact to 
100-percent contact to be small. Therefore, the usefulness of the model 
is limited to a narrow range of film thickness for any given constant load. 
In addition the location of the curves is a strong function of surface samp- 
ling interval. At this time, there is no known rational way to choose the 
sampling interval consistent with the asperity-contact phenomenon in 
elastohydrodynamie lubrication processes. 
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CONCLUSIONS 

Relations for asperity-contact time fraction as a function of nominal 
elastohydroctynamic film thickness have been presented. The calculations 
were based on a two-dimensional random surface model . Results were 
obtained for a 20-mlIlimeter-bore ball bearing with three 7. 15-milllnieter 
(9/32-in.) balls. Surface traces were obtained by using a profilometer, 
and a statistical analysis was performed in v.»hich the profile traces were 
used as statistical sample records . The results were compared with 
previously reported results based on simple one-dimensional models. 

The investigation showed that 

1. Using the two-dimensional asperity-contact model, with surface 
traces at small sample interval, the no-contact time fraction varied from 
almost full contact (90 percent) to almost no contact (1 percent) in the di- 
mensionless film thickness range 4 to 5. Full contact occurred at a film 
thickness ratio several times larger than commonly reported in the lit- 
erature. 

2. Choosing larger sample intervals in digitizing surface traces is 
the same as low- pass filtering. The effect is to lower the asperity count. 
Low-pass filtering of surface traces will shift the no-contact time frac- 
tion curves to lower values of film thickness. 

3. All the aspei'ity no-contact time fraction results may be corre- 
lated on the basis of asperity density counts, which is connected with a 
sample interval effect. 

4. A rational basis for selection of the surface sampling Interval 
consistent with the elastohydrodynamic lubrication process does not 
exist at the present time. 

5. The usefulness of the curves of contact time fraction as a means 
of determining film thickness by electrical conductance measurements is 
limited to a narrow range of film thickness ratio. 
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APPENDIX 


The probabilistic behavior of an ergodic Gaussian random surface is 
entirely defined by either the correlation function R or the power spec- 
tral density ^sd) function They are Fourier transform pairs. A 
Gaussian random process z(x,y) is one (1) that follows a Gaussian dis- 
tribution with mean p and variance and (2) for which, for all finite 
n and values of , and t , (i = l,n), the variables z(x + t 

y + T .) follow joint multivariate normal distributions. In this report 
y» * 

the correlation function is defined as 




+ T„ 




-M z(x,y)-/^^ 


(Al) 


This function measures the degree of relation between the heights of the 
random surface above two points of the reference plane that are a fixed 
distance and direction apart. 

The ergodic assumption states that the statistics of the process are 
not a function of x and y. Therefore, the expectation may be taken in 
the ensemble sense as in equation (Al; or as an average over the x-y 
plane as follows; 


rt^xpTy 

H(t , T ) - lim i- -i. / / \zpi+ T , y + T ) - pi 

y ^2 2T^ 2Ty Jt J.T L * y -i 
Ty- y 

xlz^n.y) - ^dy <Jx (A2) 


The psd is the Fourier transform of R given by 


,00 


^ y 2 K j 

(27rr 


exp E‘<Vx * Vyi ‘■’■y 


(A3) 


Hence, by means of the inverse Fourier transform the correlation func- 
tion is obtained. 
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,'SOP ,-s^oo 

/ / exp + Ty<py^ o/(<Px» 

OO ''■ -00 

The correlation function is a characterization of the surface in the 
x-y planet the psd is a characterization of the surface in the frequency 
domain, where (p and cp are frequencies. That is, a random surface 
z(x,y) may be thought of as a superposition of many surface waves of dif- 
ferent wavelength. Frequencies at which the psd is largest contribute 
more terms; frequencies at which it is smallest contribute the least. 

Spectral moments . - The interpretation of the psd as a probability 
density function leads naturally to considering the moments of that dis- 
tribution as descriptors of its shape. These moments, m.j, are defined 
here as 





x^y 


(A5) 


The psd Qjt (x,y), as defined here, is truly a probability density 
function. Thus, strictly speaking, the m^ as defined by equation (A 5) 
should not include the factor o^. The reason for including it in this way 
is so the spectral moments defined herein will correspond to the definition 
of spectral moments as presented by Longuet- Higgins and Nayak (refs. 

14 to 18). Nayak (ref. 14 ) has discussed a method of estimating the two- 
dimensional moments from a series of one-dimensional profiles obtained 
in at least three different directions. This method is explained next. 

Relation between surface and profile power spectral densities . - 
Consider a straight line through the origin in the x-y plane and at an angle 
0 to the x-axis. The height of the surface above this line is a one- 
dimensional random function of r, the distance from the origin. The 
correlation function and psd are defined by 


R0(r)=-i-E 
<>2 


IE 


(r) - n 


z (r + T) 



(A6) 


and 


Ofaitp) =~ 
^ 2n 



dr 


(A7) 
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The moments of the profile spectrum are calculated by 


(A8) 


The moments of profile psd's and the surface psd are related by the 
following equation (ref. 14): 

(A**) 

k=0 


where denotes the number of combinations of n things taken k at 
a time. 


TABLE 1. - HERTZIAN CONTACT CONDITIONS AT INNER AND OUTER 


RACES FOR THREE THRUST LOADS 


[width oi rolling track la determined by major axle width.] 


Rare 

Contact condition 

Thruat load, N (Ibf) 

80 (20) 

445 (100) 

3100 (700) 

Inner 

Maximum Hertzian atreea, GPa (kai) 
Semimajor axia, cm (in. ) 

Semiminor axia, cm (in. ) 

1.28 (185) 
0.0510(0. 0200) 
0. 0066 (0.0026) 

2.09 (303) 
0.0840 (0.0330) 
0. 0110(0. 0043) 

3.63 (527) 
0. 1500 (0. 0670) 
0.0190 (0.0076) 

Outer 

Maximum Hertzian atreaa, GRi (kai) 
Semimajor axia, cm (in. ) 

Semiminor axia, cm (in. ) 

1. 13 (164) 
0.0460(0.0180) 
0. 0086 (0. 0034) 

1.85 (269) 
0. 0740 (0.0290) 
0. 0140 (0.0055) 

fijtiHTMrjicnSfW 


TABLE n. - DIGITIZATION SAMPLE INTERVALS, NUMBER OF SAMPLE 


PCXNTS, AND TOTAL LENGTH OF SURFACE PROFILE SAMPLED 


























TABLE m. - TOTAL NUMBER OP SAMPLE PODflB 


tlBBO IN COMPUTING MOVINO AVERAGE 


Trace 

Rroaie 


Thruet load. 

N(lb) 


90 (30) 

445 (100) 

3100(700) 

1 

BaU, 0^ 

2163 

If 

Hfl 

2 

BaU, 45° 

2711 



3 

BaU. 90° 

2857 


8 143 

B 

Plat, 0° 

6623 

lot 

»7 

18 875 

B 

Flat, 18° 







B 

Plat, 36° 







mm 

Plat, 64° 







u 

Plat, 72° 







B 

Plat, 90° 







10 

Plat, 90°(r^at) 
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TABLE V. > DHflN8X>NLB88 CONI ACT FltACTlONB 


Load 

DlUMMlOlllMa 

film 

thlCkMBf, 

A 

Summits of 
htlslit A 
par squart 
eantimstar, 
KA) 

lacurakMS atwvs 
1ml A par 
squart 

rmtimstsr, 

E()0 

Oimtnslonlsss eontset fraetloaa 

N 

Ibf 

^c,ln 

T 

e,out 

^0 

Bl 

13 

4.39 

3018 

879.8 

0.717 

0.M3 

0.082 

90 


4.40 

2440 

948.9 

.573 

.879 


■a 

Ej 

4. CO 

1584 

345.7 

.384 

.429 

.801 ! 

j 

B3 

20 

4.80 

1017 

217.5 

.229 

.270 

.828 

■9 

m 

4.70 

847 

135.8 

.143 

.188 

.9298 

90 

B 

4.80 

407 

83.7 

.088 

.104 

.9727 


E 

4.90 

254 

51.2 

.094 

.084 

.9897 

■3 

20 

5.00 

157 

31.1 

.033 

.039 

.99822 

■3 


5.10 

98 

18.7 

.020 

.023 

.99883 

na 

100 

4.55 

1242 

287.8 

0.771 

0.888 

0.037 i 

E9 
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4.80 

994 
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.810 
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.380 
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.289 

.284 

.826 
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3100 700 


9.30 


30.3 

16.2 


82.8 

90.8 

30.7 
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KiS 

0.727 
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BSI 

.444 

.498 

IkI 
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.302 
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.9893 

.057 
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Figure t - Croti section of e cow p e s He surhee rou^neis precau. 
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